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Searches for low-energy tetrahedral polymorphs of carbon and silicon have been performed using 
density functional theory computations and the ab initio random structure searching (AIRSS) ap¬ 
proach. Several of the hypothetical phases obtained in our searches have enthalpies that are lower or 
comparable to those of other polymorphs of group 14 elements that have either been experimentally 
synthesized or recently proposed as the structure of unknown phases obtained in experiments, and 
should thus be considered as particularly interesting candidates. A structure of Pbam symmetry 
with 24 atoms in the unit cell was found to be a low energy, low-density metastable polymorph in 
carbon, silicon, and germanium. In silicon, Pbam is found to have a direct band gap at the zone 
center with an estimated value of 1.4 eV, which suggests applications as a photovoltaic material. We 
have also found a low-energy chiral framework structure of P4i2i2 symmetry with 20 atoms per cell 
containing fivefold spirals of atoms, whose projected topology is that of the so-called Cairo-type two- 
dimensional pentagonal tiling. We suggest that P4i2i2 is a likely candidate for the structure of the 
unknown phase XIII of silicon. We discuss Pbam and P4i2i2 in detail, contrasting their energetics 
and structures with those of other group 14 elements, particularly the recently proposed P42/ncm 
structure, for which we also provide a detailed interpretation as a network of tilted diamond-like 
tetrahedra. 

PACS numbers: 71.15.Nc,62.50.-p,61.50.-f,61.66.-f 

I. INTRODUCTION 

The group 14 elements carbon (C), silicon (Si), and 
germanium (Ge) have attracted much interest and have 
been extensively studied. These elements have an s 2 p 2 
valence electronic configuration, which leads to common 
chemical features, but also significant differences. Pure 
C is found on the Earth mainly in the graphite and cubic 
diamond (cd) forms, which exhibit some of the strongest 
bonds known in Nature. Pure Si and Ge also adopt the 
diamond structure under ambient conditions, with ideal 
tetrahedral coordination, and they are, of course, semi¬ 
conductors of great importance in the electronics indus¬ 
try. Under applied pressure of around 11 GPa both Si 
and Ge transform into the /3-Sn-type structure^ which 
has 4+2 coordination and is metallic. Other stable and 
well researched phases of these elements exist at higher 
pressures— *2 

Besides their thermodynamically stable forms, a num¬ 
ber of zero-pressure metastable polymorphs have been 
observed, with both higher and lower densities than 
the corresponding diamond phases. The lonsdaleite or 
hexagonal diamond form of carbon, which has a density 
similar to cubic diamond, was first observed in mete¬ 
orite craters^ and later synthesized from graphite under 
high pressure and temperature^ and also by chemical 


vapor deposition and other chemical methods. Hexag¬ 
onal diamond has also been reported in Si and Ge af¬ 
ter decompression and heating from the high-pressure 
phases or after indentation— - — Among the high-density 
polymorphs, the cubic bc8 phase of Si and tetragonal 
stl2 phase of Ge were first obtained about fifty years ago 
as zero-pressure tetrahedral metastable forms recovered 
upon decompression from the respective high-pressure f3- 
Sn-type phases— The r8 structure (a rhombohedral 
distortion of the bc8 structure) was later discovered ex¬ 
perimentally in high-pressure experiments on quenched 
Si™— and has recently been identified in nanoindented Ge 
samples — - Although the bc8 structure has also been re¬ 
ported in Ge ) 12 ’ 13 the stl2 structure has not so far been 
observed in Si—- Low density open framework clathrate- 
type structures of Si and Ge, with several potential appli¬ 
cations, can also be produced at ambient pressure using 
chemical synthesis methods and have been the object of 
intense research in recent years—There are also sev¬ 
eral reports of other Si and Ge phases obtained upon 
pressure release in diamond anvil cells or by nanoindenta¬ 
tion and whose structure has not yet been experimentally 
resolvedi 6 ! 7 : 19- — 

Carbon, with the versatility that stems from its unique 
ability to form sp , sp 2 , and sp 3 hybrid bonds, pro¬ 
vides most examples of allotropes with a wide range 
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of structures and properties. These allotropes include, 
besides those already mentioned, fullcrenes, nanotubes, 
and graphene and its close relatives.— New structural 
forms have been added to this list in recent times. 
A low-temperature quenchable transparent and super- 
hard metastable crystal phase of carbon, for which sev¬ 
eral structures have been proposed, can be obtained 
by cold compression of graphite*- 24 whereas the applica¬ 
tion of pressure under different controlled conditions to 
fullerenes and carbon nanotubes leads to the formation 
of some new and elusive forms——' These findings have 
excited the imagination and attracted the attention of 
researchers to this field and it is expected that new ex¬ 
perimental observations will occur in the coming years. 

Along with experimental investigation of observed al¬ 
lotropes, great effort has been expended in theoretical 
searches for new tetrahedral phases whose predicted me¬ 
chanical, structural and optical properties may have po¬ 
tential in advanced technological applications. Such the¬ 
oretical experiments may result in a fruitful interplay 
between theory and practice. There is a need for com¬ 
puter simulations that help in characterizing new phases 
obtained in the laboratory precisely, as such phases are 
difficult to study in high-pressure experiments, and sim¬ 
ulations can significantly help in identifying the correct 
structure, or in discarding others. 

A number of novel hypothetical phases have been pro¬ 
posed from these theoretical studies. One driving force in 
the search for new carbon phases has been the quest for 
superhard materials - novel materials with hardness ri¬ 
valling or even exceeding that of carbon diamond, which 
could be of technological importance— r— For silicon, 
the main material used so far in the fabrication of solar 
cells, there is interest in new phases with electronic band 
structures and optical properties better suited than the 
diamond phase for photovoltaic applications! 33 ' 34 There 
is also a purely fundamental interest in researching the 
possibilities of structural and functional diversity of this 
important class of elemental materials. 

Here we describe a number of hypothetical energet¬ 
ically favourable polymorphs of C, Si, and Ge, in¬ 
cluding two particular low-energy polymorphs of space 
group symmetries Pbam and P4i2i2. These structures 
were obtained in density-functional-theory (DFT) based 
searches. They are somewhat less dense than the corre¬ 
sponding diamond structures and have very similar den¬ 
sities to a low-energy P42/ncm structure that was re¬ 
cently found in DFT structure searches^ and has been 
considered as a candidate for the observed metastable 
phase XIII of Si£— —~— The Pbam structure has a some¬ 
what lower energy in C, Si, and Ge than the P4i2i2 and 
P4 2 /ncm structures. It may be hoped that some of these 
low-energy polymorphs can be synthesized experimen¬ 
tally. 

The rest of the paper is organized as follows. In Sec. HD 
we give details of our ab initio calculations and the struc¬ 
tural searches performed. In Sec. m we show our results 
for the materials considered, including detailed descrip¬ 


tions of several selected structures. Finally, in Sec. [V] we 
present a summary of our most important conclusions. 
The plentiful data generated during the study has made 
it necessary to place a portion of them as Supplemental 
Material to the present paper— 


II. AB INITIO CALCULATIONS AND 
STRUCTURE SEARCHING 

We performed computational searches for low-energy 
structures of C and Si using first-principles density- 
functional-theory (DFT) methods and the ab initio ran¬ 
dom structure searching (AIRSS) approach— ^ We did 
not perform searches for Ge, instead we took the low- 
energy Si structures and re-relaxed them for Ge, as low- 
energy polymorphs of Si are also expected to be low- 
energy polymorphs of Ge. In the AIRSS approach an 
ensemble of randomly chosen initial structures are re¬ 
laxed to a minima of the energy. AIRSS has been suc¬ 
cessfully employed in finding low-energy structures in 
many systems, including group 14 elements and their 
compounds— r— We performed searches for structures 
with up to 24 atoms per unit cell. 

The CASTEB^ DFT code and the Perdew-Burke- 
Ernzerhof (PBE)^ Generalized Gradient Approximation 
(GGA) density functional were employed for the searches. 
We used ultrasoft pseudopotentials^I with the outermost 
four valence electrons treated explicitly and default basis 
set energy cutoffs, further relaxing the structures of in¬ 
terest at a higher level of accuracy, using a k-point grid 
of spacing 27 tx0.03 A -1 for the Brillouin zone integra¬ 
tions. The details of this procedure are similar to those 
reported in previous studies—^ 

The structures were subsequently re-relaxed using the 
PBE functional and the Vienna ab initio simulation pack¬ 
age (vasp)j 48 ' 49 For Ge we tested pseudopotentials in¬ 
cluding the four valence electrons explicitly and pseu¬ 
dopotentials in which the 3d electrons of Ge were also in¬ 
cluded explicitly, but the effects of including the 3d elec¬ 
trons were very small, and we chose to perform phonon 
calculations with the d electrons treated as core states. 
Calculations of the energies as a function of volume, 
phonon spectra, and band structures were performed us¬ 
ing the projector augmented-wave (PAW) methodi 50 ' 51 
The phonon calculations were performed with the PHON 
code— The plane wave cutoff energies used for the VASP 
calculations were 520 eV (C), 320 eV (Si), and 230 eV 
(Ge, 375 eV for the calculations with the 3d electrons). 
Dense k-point grids were used for the Brillouin zone in¬ 
tegrations (e.g., an 8x8x8 grid was used for the Pbam 
structure, with 24 atoms per cell). The structural relax¬ 
ations were deemed to be converged when all of the forces 
were less than 5 meV A -1 and the anisotropy of the stress 
tensor was less than 0.1 GPa. The results of the CASTEP 
and VASP calculations were in excellent agreement with 
one another, which is a useful crosscheck. We also per¬ 
formed VASP calculations using the local spin-density ap- 
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proximation (LSDA)^ and PBEsol density functionals^ 
which gave qualitatively very similar results to the PBE 
calculations and allowed us a full comparison among the 
three functionals. 

We calculated the total energy E as a function of vol¬ 
ume V for each phase and fitted the data using a fourth- 
order Birch-Murnaghan equation of stated (EoS) from 
which the pressure, p, and the enthalpy, H = E + pV, 
were obtained. 


III. 

RESULTS 

A. 

Energetics 


Enthalpy-pressure relations for a selection of different 
structures in C, Si, and Ge are shown in Figs. [TJ [2j and 
[U Calculated equilibrium zero-pressure volumes Vo, bulk 
moduli B 0 , and pressure derivatives of the bulk modulus 
B' 0 , obtained from the EoS fitting, are collected in Tables 
UTlIIIIl and IIVI for C, Si, and Ge, respectively. The differ¬ 
ences in energy from the corresponding diamond struc¬ 
tures at zero pressure are also reported in Tables ITIllIVI 
For a given structure, these energy differences generally 
decrease with increasing atomic number. (See also Sup¬ 
plemental Material— i) 

In general, we find good agreement between the cal¬ 
culated and experimental values of structural parame¬ 
ters for the observed phases. For example, our values 
of the volume, bulk modulus and pressure derivative of 
the bulk modulus of C-cd are in good agreement with 
the experimental values of Vo = 5.6738(13) A 3 (room 
temperature)^ B 0 = 444(4) GPa£! and B' 0 = 3.65(5), 
data from similar DFT calculations ) 581 ’ 59 and values from 
accurate quantum Monte Carlo calculations , 60 ' 61 We 
have adopted the experimental data of Ref. H3 for B 0 
and Bq of C-cd, with the revised pressure scale of Ref. 
l58l . Likewise, we obtain good agreement with the exper¬ 
imental data for Si-cd and Ge-cd. 


B. Overview of the findings 

It is normally possible to identify families of related 
crystalline structures. For example, the cd and hd struc¬ 
tures are members of an infinite family of low-energy 
polytypic structures consisting of different stackings of 
layers. Structures consisting of (1) different stackings of 
layers; (2) regions of the most stable cd phase and interfa¬ 
cial regions; (3) periodic arrays of point defects in the cd 
structure etc., can be dreamt up with energies per atom 
arbitrarily close to that of the corresponding cd struc¬ 
ture. It is not our intention to investigate such families 
of tetrahedrally-bonded group 14 structures, but instead 
we seek structures which differ in substantial ways from 
known structures. 

Our searches produced many structures, including the 
observed equilibrium cubic diamond (cd) structure of C 


and Si, the hexagonal diamond structure (hd)^ and 
dense structures adopted by Si under high pressures, 
such as the /3-Sn structure. As well as finding experi¬ 
mentally observed tetrahedrally-bonded phases we found 
many structures that have been reported in recent DFT 
searches of group 14 elements and considered others ex¬ 
plicitly for the sake of compariso m 29 i 35i42 ^~ 6 - (Though 
we will not report on all of them, some further data can 
be found in the Supplemental Material.—) We refer to 
the new structures by the names of their space groups 
and use common abbreviations for the well established 
structures, although some structures will be occasionally 
also referred to by other names with which they have ap¬ 
peared in the recent literature on this subject. We focus 
here on tetrahedrally bonded low-energy structures. 


1. High-density tetrahedral polymorphs 

In our searches we found all the dense tetrahedral 
polymorphs that have been observed in Si and/or Ge: 
the bc8 and r8 structures obtained upon decompres¬ 
sion from the corresponding high-pressure /3-Sn phase or 
nanoindentation^iSr— iSI and the stl2 structure which 
has only been observed upon decompression in Ge4 i 13 i 14 
We also found a body-centred tetragonal structure with 
space group /4i/a and 8 atoms per primitive cell, whose 
enthalpy is remarkably close to that of r8/bc8, and which 
plays a role in the decompression kinetics from the high 
pressure phases—2. For Si within the PBE, for example, 
I4i/a is only 5 meV per atom higher in enthalpy than 
bc8/r8 at zero pressure (and about the same enthalpy 
as stl2, within the accuracy of our calculations), see Ta¬ 
ble IIIII which makes I4i/a an energetically competitive 
high-density polymorph. 

Also belonging with the group of high-density tetrahe¬ 
dral polymorphs, an hexagonal structure of P 6422 sym¬ 
metry labelled cintet is interesting as it is the elemental 
analog of the tetrahedral variant of the binary cinnabar 
structure (also called pseudo-cinnabar structure) ob¬ 
served at high pressures in both ZnTe and GaAs—At 
zero pressure, the cintet structure in Si and Ge is very 
close in energy to the Ibam structure with 6 atoms per 
primitive cell recently proposed by Malone and Cohen^ 
They are both, however, higher in energy than the other 
observed dense allotropes bc8/r8/stl2, which may im¬ 
pede their direct synthesis as metastable phases on pres¬ 
sure release from the high-pressure /3-Sn phase. At higher 
pressures Ibam is lower in enthalpy than cintet. In gen¬ 
eral, the high-density allotropes have equilibrium ener¬ 
gies significantly higher than those of the low-density 
forms that we discuss next. However, their reduced vol¬ 
umes make them competitive at high pressures, and in¬ 
deed it is under such conditions that some of them have 
been synthesised. Among the phases considered, the cin¬ 
tet phase has a bulk modulus second only to the diamond 
forms of C, and somewhat higher to their respective dia¬ 
mond forms in Si and Ge. 
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2. Low-energy and lower-density tetrahedral polymorphs 

One interesting result of our search is that we have 
found a C-centred orthorhombic structure of Cmca sym¬ 
metry containing 16 atoms in the primitive unit cell 
whose E(V) curve is almost degenerate, in each of the 
group 14 elements studied, with that of the Cmmm 
structure (with 8 atoms per primitive cell) previously 
found by Zhao et air 0 in carbon and named Cco-Cs 
by these authors. This Cmmm structure has been pro¬ 
posed as that of a superhard carbon allotrope recovered 
from room-temperature compression of carbon nanotube 
(CNT) bundles^ and it is indeed structurally related to 
the (2,2) CNT.— The Cmmm carbon phase of Zhao et 
al. was also readily found with AIRSS. 

From our calculations, both the Cmmm and Cmca 
structures are very close in energy to the Pmmn structure 
of the so-called P carbon phase, with 16 atoms per cell, 
alternatively proposed by Niu et as the structure of 
the cold-compressed CNT obtained in experiment in Ref. 
l27l . Both the new Cmca allotrope and the Cmmm and 
Pmmn allotropes share similar structural motifs, namely 
variations on arrangements of the same four-, six-, and 
eightfold ring patterns, which is the origin of their simi¬ 
lar energetics. They also contain many sixfold diamond- 
like rings which may partially account for their relatively 
low energy and superhardness in carbon. Possibly, other 
energetically close structures could be conceived exhibit¬ 
ing variations of these motifs. This suggests that cold 
compressed CNT bundles could adopt a mix of different 
but structurally and energetically very close forms, as 
described by the family Cmmm/Pmmn/Cmca. 

Our search also produced the primitive tetragonal 
structure with 12 atoms per cell (T12) and space group 
P^/ncm that has been recently reported by Zhao et 
al£ 5 In the present work we provide a detailed descrip¬ 
tion of this structure and its relation to diamond (see Sec. 
mm , interpreting P42/ncm as a network of diamond- 
like tetrahedra made up of slabs of almost regular corner¬ 
sharing tetrahedra alternatively rotated left and right, so 
that fivefold connectivity appears within each slab and 
adjacent slabs are bonded together by non-shared cor¬ 
ners. We believe that this remarkable structure is the 
simplest that can be built by the stacking of slabs of 
tilted tetrahedra while preserving to a large degree a 
highly regular fourfold coordination for the sites. The 
P42/ncm structure has been proposed^ as a candidate 
for both the Si-XIII phase (which is observed to coexist 
along with the Si-cd, Si-hd, Si-r8, and Si-bc8 phases in 
nanoindentation experiments^^ - — ) and for an experi¬ 
mentally synthesized metastable Ge phased and for this 
reason we compare our data for other candidate phases 
with those of P42/ncmJ 1 

We have also obtained two interesting and novel low- 
energy structures, Pbam (with 24 atoms per cell) and 
P4i2i2 (20 atoms per cell), with orthorhombic and 
tetragonal symmetry, respectively, which are energeti¬ 
cally competitive polymorphs and whose detailed de¬ 


scription will be given below. The orthorhombic Pbam 
structure and the tetragonal P4i2i2 structure are slightly 
expanded in volume with respect to the corresponding 
diamond structures, with zero-pressure equilibrium vol¬ 
umes about 2-3 % larger than those of the correspond¬ 
ing cd structures, see Tables HTHlYI The density of the 
P 6522 chiral framework structure reported previously^ 
is intermediate between those of Pbam/P4i2i2/P42/ncm 
and the type-I and II clathrates. The calculated phonon 
dispersion curves of these phases show that they are dy¬ 
namically stable. 

Note that the present Pbam structure with 24 atoms 
per cell is not the same as the Pbam structure with 16 
atoms per cell proposed by Niu et al. for R carbon^ 
along with other orthorhombic structures. It is in par¬ 
ticular significantly lower in energy. The present Pbam 
structure is energetically highly favorable in C, Si, and 
Ge. For carbon, Pbam is essentially degenerate in en¬ 
ergy with the most favourable of the low-symmetry defect 
structures proposed recently by Botti et al ., with mon¬ 
oclinic symmetry and 20, 22, and 24 atoms per cell.' 2 
whereas in our study for silicon and germanium we have 
considered ex professo these hypothetical carbon struc¬ 
tures and have found them to lie within 4 meV of Pbam, 
which is close to the accuracy of the calculations— Apart 
from these structures, among those considered here, only 
the diamond forms and the type-II clathrate of germa¬ 
nium (Ge-clat34) are significantly lower in enthalpy at 
zero pressure than the corresponding Pbam structures, 
with the type-I germanium clathrate (Ge-clat46) having 
about the same zero-pressure entahlpy as Ge-Pbam (and 
a rather smaller density). 

Low density clathrate structures of types I and II 
(here labelled clat46 and clat34, respectively) have been 
synthesized in Si, Ge, and Sn^ - — although not in 
C. Group 14 clathrates have normally been synthe¬ 
sized by incorporating guest atoms, although a guest- 
free Ge clathrate has been formed.— Our calculated re¬ 
sults for the clathrates are in good agreement with earlier 
work. 74- ' 8 The high-density polymorphs bc8, r8 and stl2 
of Si and Ge are all considerably less stable than Pbam 
at zero pressure. 

The P4i2i2 form is a chiral framework structure con¬ 
taining fivefold spirals of atoms, that is also energeti¬ 
cally competitive, particularly in Si, for which among 
the new polymorphs only the Pbam structure is lower in 
energy. In C, its enthalpy is very close to those of the 
Cmmm/Cmca phases already mentioned in connection 
with cold compressed nanotubes. The P42/ncm struc¬ 
ture is less stable than P4i2i2 in Si, but more stable 
than P4i2i2 in C whereas for Ge both structures have ap¬ 
proximately the same enthalpy. The two pure clathrate 
structures considered are also less stable than P4i2i2 in 
Si. 

As has been observed previously, the calculated energy 
differences within the PBE are systematically larger than 
those from the PBEsol functional, and these are larger 
in turn than those from the LDA. One effect of this is 






5 


that the coexistence pressures with the diamond phases 
for the high pressure phases (such as /3-Sn) or the high- 
density tetrahedral metastable phases (such as bc8) are 
larger with PBE than PBEsol or the LDA. The differ¬ 
ences in energy between the PBE and the LDA or PBEsol 
are, however, quite small for the low-density tetrahedral 
polymorphs with equilibrium volumes close to the dia¬ 
mond phase such as Pbam, P42/ncm, or P4i2i2, whereas 
for the high-density tetrahedral polymorphs such as bc8, 
r8, or stl2 there is a significant difference between the 
results obtained with different exchange-correlation ap¬ 
proximations, see Tables and Supplemental Mate¬ 
rial. This reflects the similarity in bonding in diamond 
and in phases such as Pbam, P42/ncm, or P4i2i2. The 
bulk moduli of these and the other low-density struc¬ 
tures obtained in the searches tend to be slightly below 
the bulk moduli of the corresponding diamond forms. In 
the case of carbon, for example, these are in the range 
~400-415 GPa (cf. the calculated value for C-cd of 435 
GPa). 



Pressure (GPa) 

FIG. 1. (Color online) Enthalpy-pressure relations for various 
structural phases of carbon calculated using the PBE. The 
enthalpies are given with respect to the zero-pressure C-cd 
phase. 


C. Description of selected low-energy structures 

The three phases on which we focus next, P42/ncm, 
Pbam, and P4i2i2, are structurally very different among 
themselves and with respect to the diamond structure. 
Topologically, all the atoms in these structures have a 
tetrahedral coordination of nearest neighbours, in some 
cases considerably distorted from a regular coordination, 
presenting five-, six-, and/or sevenfold rings of atoms. 
Each of these novel structures has similar characteristics 



FIG. 2. (Color online) Enthalpy-pressure relations for various 
structural phases of silicon calculated using the PBE. The 
enthalpies are given with respect to the zero-pressure Si-cd 
phase. 




Pressure (GPa) 


FIG. 3. (Color online) Enthalpy-pressure relations for various 
structural phases of germanium calculated using the PBE. 
The enthalpies are given with respect to the zero-pressure 
Ge-cd phase. 


in all four materials studied. Hereafter we will discuss 
them in detail providing, for the sake of illustration, cal¬ 
culated numerical data (values of distances, angles etc) 
which correspond to Si at zero pressure, unless otherwise 
stated. In Table [T] we summarize the relevant crystal¬ 
lographic structural data for the P4 2 /ncm, Pbam, and 
P4i2i2 structures of Si, at zero pressure, calculated us¬ 
ing the PBE functional. Diffraction patterns of the struc- 
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TABLE I. Crystallographic data for the I4i/a, Cmca, 
P42/ncm, Pbam, and P4i2i2 structures of Si at zero pressure, 
calculated using the PBE: space group (SG), lattice param¬ 
eters, and atomic sites (Wyckoff positions). Those internal 
coordinates fixed by symmetry are given as integer fractions, 
to aid recognition. 


I4i/a 

SG Mi/a, No. 88 
a=6.676 A, c=6.514 A 

Sil 16f (-0.0962, -0.0994, 0.1170) 

Cmca 

SG Cmca, No.64 


a=6.480 A, 6= 

15.496 A, c=6.693 A 


Sil 

16g 

(0.3161, 0.0673, 0.4139) 


Si2 

16g 

(0.3154, 0.1887, -0.0849) 

P4 2 /ncm 

SG P^/ncm, No. 138 


a=5.221 A, c= 

9.295 A 


Sil 

4b 

(3/4, 1/4, 1/4) 


Si2 

8i 

(0.0865, 0.0865, 0.3924) 

Pbam 

SG Pbam, No. 

55 


a=11.774 A, b- 

=11.081 A, c=3.863 A 


Sil 

4h 

(0.6158, 0.2400, 1/2) 


Si2 

4h 

(0.4625, 0.1005, 1/2) 


Si3 

4h 

(0.9081, 0.9569, 1/2) 


Si4 

4g 

(0.7705, 0.7075, 0) 


Si5 

4g 

(0.8040, 0.0124, 0) 


Si6 

4g 

(0.1501, 0.6240, 0) 

P4i2i2 

SG P4i2i2, No. 92 


a=8.903 A, c= 

5.261 A 


Sil 

8b 

(0.4312, 0.1390, 0.7001) 


Si2 

8b 

(0.4752, 0.6504, 0.4205) 


Si3 

4a 

(0.2011, 0.2011, 1/2) 


tures simulated using the fullprof software- 9 , are shown 
in Fig. 0 Further details of these structures are given in 
the Supplemental Material^ 

1. The Pfa/ncm structure as a network of tilted tetrahedra 

P42/ncm is a remarkably simple structure that can 
be understood as a stacking of (001) slabs of diamond¬ 
like corner-sharing tetrahedra, with adjacent tetrahedra 
within each slab alternately tilted clockwise and anti¬ 
clockwise by approximately 20°, see Fig. [5] The cen¬ 
ters of the tetrahedra correspond to the Sil sites of the 
P4 2 /ncm structure (red sites in Fig. [5]) and their corners 
correspond to the Si2 sites (blue). There are only three 
free internal parameters in P4 2 /ncm, all them related to 
the Si2 corner sites, whereas the Sil center sites form a 
simple tetragonal lattice. Within each of these tetrahe- 



20 (deg) 


FIG. 4. (Color online) Simulated powder XRD patterns of 
the P42/ncm, P4i2i2, and Pbam phases of Si at zero pres¬ 
sure, calculated within the PBE approximation. The x-ray 
wavelength used was 0.44864 A. 


dral slabs, the rotation of adjacent tetrahedra in opposite 
sense brings their corners together, with the formation of 
new bonds between them (blue rods in Fig. 03), above 
and below the plane of the centers and orientated along 
(110) directions. When two of these tetrahedral slabs, 
oppositely orientated (Fig. 03 and c), are stacked along 
the z axis, new inter-slab bonds are established between 
the corners of each tetrahedron and the corners of the op¬ 
posite tetrahedra immediately above and below it (green 
rods in Fig. 01). Note that the centers of the tetrahedra 
in different slabs are strictly aligned along the stacking 
direction, whereas the orientation of the stacked tetra¬ 
hedra are alternately tilted clockwise and anti-clockwise. 
As a result of intra- and inter-slab bonding of tetrahe¬ 
dral corners, each Si2 corner site is linked to two other 
Si2 corner sites (as well as to two Sil center sites), with 
the formation of staggered chains of Si2 corner sites run¬ 
ning along (110) directions. 

It is quite remarkable that while the tilting and rebond¬ 
ing reduces the symmetry of the P4 2 /ncm tetrahedral ar¬ 
rangement with respect diamond, so that in P4 2 /ncm the 
tetrahedral units become tetragonal disphenoids, these 
units are actually very close to regular tetrahedra. The 
Sil center sites have an almost ideal tetrahedral Si2 en- 
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vironment, with all four bonds of equal length and near¬ 
ideal tetrahedral angles. The tetrahedral units perform 
largely as rigid units in the P 42 /ncm structure. 

The tilting and subsequent intra-slab corner-bonding 
leads to the appearance of fivefold rings made up of two 
Sil tetrahedra centers and three Si2 corners, see Fig. 
G3>- Each of these fivefold rings is a slightly buckled 
and symmetrical pentagon with four out of its five sides 
having the same length, and with three different angles, 
which deviate only moderately from the angle of 108° of 
a regular pentagon. The stacking of opposite tetrahedral 
slabs and accompanying corner-bonding further results 
in twisted sixfold rings ( 2 xSil, 4xSi2) with alternating 
orientation along the z axis, noted by Zhao et al and 
of sevenfold rings ( 2 xSil, 5xSi2) in a chair-like configu¬ 
ration, also alternating in orientation along the stacking 
direction. The twisted sixfold rings are quite unlike those 
in the diamond structure. There is only one type of five-, 
six-, and sevenfold ring. 

There are only three different bond lengths in 
P 42 /ncm: one (di) between Sil and Si2 sites (that is, 
from the center to each corner of the quasi-regular tetra¬ 
hedra), another between neighboring Si2 sites in the same 
slab (cfe), and a third between Si2 sites in adjacent slabs 
(dz). These three different bond lengths are represented, 
respectively, by red, blue, and green rods in Fig. [5] The 
tetrahedral distance d\ is very similar to that of the dia¬ 
mond phase (2.37 A) and shows a similar compressibility 
under applied pressure, which is in line with the fact 
that the tetrahedral units of P 42 /ncm remain very much 
diamond-like in nature. The inter-slab distance d$ is also 
similar though it has a slightly larger compressibility, 
while the intra-slab distance g ?2 is larger than both d\ 
and c ?3 and it also shows a larger compressibility. There 
are also six different bond angles^ with only one angle 
deviating significantly (122.5°) from the ideal tetrahe¬ 
dral value (109.47°). This large angle corresponds to the 
chains of bonds between Si2 tetrahedra corners running 
along (110) directions. The bond angles in P 42 /ncm vary 
very little when the structure is compressed. The large 
Si2-Si2-Si2 angle increases by about 1.5° in the range 
from -5 GPa to 10 GPa, while the rest of the angles 
change by tenths of a degree. The tilting of the tetra¬ 
hedra is locked by intra- and inter-slab bonding and the 
aspect of the P4 2 /ncm structure would not change much 
even over a substantial pressure range. 


2. The Pbam structure 

With six inequivalent sites and twelve free internal 
parameters, Pbam has a larger variation of tetrahedral 
environments than P 42 /ncm, and yet many elements of 
Pbam remain remarkably similar to those of diamond 
with which a useful comparison can be established. In 
Fig. [G] we represent the Pbam structure in a way that 
allows comparison to the diamond structure, when the 
latter is viewed along its (110) direction. In analogy to 



FIG. 5. Construction of P 42 /ncm as a network of tilted 
diamond-like tetrahedra. Panel (a) shows a single slab of 
corner-sharing regular tetrahedra of the diamond structure 
viewed along the (001) direction. The projected tetrahedra 
and conventional cubic cell are shown in black. The centers 
of the tetrahedra are shown in red and their corners in blue, 
with the different sizes of the atoms and varying thickness of 
the bonds used to give a sense of perspective. Panels (b) and 
(c) show two adjacent slabs (lower and upper, respectively) 
of tilted tetrahedra that alternate stacked along the 2 axis 
to give the P 42 /ncm structure, as shown in (d)j£H Note the 
opposite orientation of the tetrahedral bonds and the tilting 
of the tetrahedra in the two adjacent slabs. Intra-slab bonds 
(parallel to the plane of projection), resulting from the tilting 
of the tetrahedral units within each slab, are represented by 
blue rods, while inter-slab bonds, arising from the stacking of 
the slabs, are represented by green rods. Intra-slab bonding 
leads to characteristic fivefold rings, one of which is depicted 
in blue in (b), while inter-slab bonding leads to buckled sixfold 
rings, two of which are shown in profile in (d). 


diamond, we see stacks of zigzag chains in Pbam, two per 
cell, linking the three 4g-type sites Si4-6 (red sites in Fig. 
0 ). These distorted diamond-like chains are contained in 
( 001 ) projection planes and run along the ( 010 ) direction. 
The three 4h-type sites Sil-3 (blue sites) form in turn two 
staggered planar chains along the same direction but at 
different (001) planes from the 4g-chains and intercalated 
between them. Also note the relative displacements and 
different orientations of consecutive 4g and 4h chains. 
Adjacent 4g and 4h chains are connected by three differ¬ 
ent types of bond which form zigzag diamond-like chains 
running perpendicular to the projection plane, each link¬ 
ing one type of 4g site to one type of 4h site (Sil-Si 6 , Si2- 
Si4, and Si3-Si5), so each 4g site (respectively, 4h site) 
is linked to two 4g (4h) sites on the same (001) plane 
and to two 4h (4g) sites of the same type, above and be¬ 
low the plane. The orientations of these three different 

























connecting 4g-4h chains of bonds deviate greatly from 
the regular orientation found in the diamond structure 
(which would correspond to a (010) plane in Pbam) but 
the chains themselves are very little distorted compared 
with equivalent chains in diamond. 

The main topological difference between diamond and 
Pbam arises from the intra-chain connectivity within the 
4h slabs in which drastic changes lead to the formation 
of the fivefold and sevenfold rings shown in Fig. [ 6 ] There 
are in fact no other major topological changes either 
within the 4g slabs or within the interconnecting 4g-4h 
chains. A simple interpretation of Pbam is that it con¬ 
sists in a certain pattern of shearing of atomic slabs in 
diamond such that breaking and formation of intra-slab 
bonds occurs only within alternate slabs (the stacks of 4h 
chains in Pbam) whereas the connecting chains of inter¬ 
slab bonds are not affected beyond a small deformation 
of bond lengths and angles. All in all, two bonds are 
broken and reformed per 4h chain and cell, which drasti¬ 
cally transforms the 4h intra-slab connectivity giving rise 
to fivefold and sevenfold rings perpendicular to the shear¬ 
ing direction, though some sixfold diamond-like rings are 
also preserved, see Fig. [ 6 ] Whereas the 4h slabs undergo 
an important reconstruction, the topology of the 4g slabs 
remains essentially unchanged (though deformed) and 
the 4g-4h chains of connecting bonds are mostly unal¬ 
tered (though rotated from their initial orientation by 
the intra-4h rebonding). These changes are accompanied 
by a contraction in the y direction and elongation in the 
x direction with respect to diamond. 

There is one type each of fivefold, sixfold, and seven¬ 
fold rings in the xy projection of Pbam shown in Fig. O 
which appear in four different orientations 42. Apart from 
the sixfold rings shown in Fig. [ 6 j the Pbam structure 
preserves a large number of other sixfold chair-like rings 
of the diamond structure perpendicular to the xy plane 
which are shown in profile in this figure. (There are also 
eightfold rings perpendicular to the xy plane formed by 
the rebonding of the 4h slabs). 

Due to the large number of free parameters of the 
Pbam structure, there are ten different nearest neigh¬ 
bor distances and twenty one different bond angles. Dis¬ 
tances vary between 2.33 and 2.42 A (cf. the calculated 
value for the diamond structure of 2.37 A) with most 
distances within 0.4 % of the diamond value and a mean 
distance of 2.37 A. Bonds within the 4g slabs are some¬ 
what shorter than diamond bonds whereas those within 
the 4h slabs are somewhat longer, with the longest bond 
occurring between the Si2 and Si3 sites that form the 
pentagonal rings. Bond compressibility is also similar to 
that in diamond, although the two Sil-Si2 and Si3-Si3 
bonds within the 4h slabs have somewhat larger com¬ 
pressibility and one Si5-Si6 bond within the 4g slabs has 
a smaller compressibility. 

There is also a considerable angular distortion with re¬ 
spect to the ideal tetrahedral angle. Not surprisingly, 
the largest angle (123.9°) occurs within the heptagonal 
rings (two other different angles in the sevenfold rings 



FIG. 6. (Color online) A slab of the Pbam structure viewed 
along the 2 axis, with unit cells shown in black and a modicum 
of perspective so as to reveal the connectivity among sites. 
The six inequivalent sites are labelled and shown in different 
colors: the 4g-type sites are depicted in three different hues 
of red, while blue is used for the three 4h-type sites. The 
fivefold and sevenfold rings of the structure (as well as sixfold 
diamond-like rings) are all visible in this projection. Note 
that the red 4g sites and the blue 4h sites are each in different 
projection planes a distance c/2 apart. 


are close by, at 120 . 6 °; cf. the angle of a regular hep¬ 
tagon of 128.57°). The minimum angle (92.6°) occurs 
at the vertex of the fivefold rings (where the rest of the 
angles vary between 101.4° and 102.7°). In spite of these 
distortions it is quite remarkable that the Sil sites at the 
junction of two adjacent sixfold rings have a tetrahedral 
environment of neighbors very similar to that of diamond 
sites. In fact, the sixfold rings themselves are very similar 
to those present in ideal diamond, with angles deviating 
slightly from the ideal tetrahedral angle (from 107.7° to 
110.6°) and a mean angle of 109.3°. The angular devi¬ 
ations with respect to diamond are also minimal along 
the three different zigzag chains running perpendicular 
to the xy plane (108.6°, 109.2°, and 110.1°) where the 
bond lengths are also very similar to diamond. As in 
P 42 /ncm, the bond angles remain very stable under com¬ 
pression with a maximum variation of 1.8° from -5 GPa 
to 10 GPa, but in fact most angles vary only by tenths 
of a degree. 


3. The P4i 2\ 2 structure: A new chiral framework with 
tetragonal symmetry 

Whereas Pbam and P 42 /ncm share some common ele¬ 
ments with the diamond structure, P4i2i2 does not show 
any obvious similarity with diamond, beyond being also 
a tetrahedral network. With three different types of site, 
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six different bond lengths and sixteen different bond an¬ 
gles, such a network is quite unlike the regular tetrahedral 
network found in diamond. Instead, P4i2i2 has an ob¬ 
vious similarity to the chiral framework structure (CFS) 
of P 6522 symmetry previously discovered by Pickard and 
Needs . 42 i 83 

In Fig. [7] we show the P4i2i2 structure viewed from 
the direction of its main rotation axis (cf . the Supple¬ 
mental Material^ and Fig. 1 of Ref. |42| for a similar 
plot of the CFS-P 6522 structure). The two-dimensional 
pattern of this projection corresponds to the so called 
Cairo pentagonal tiling, which is reported to appear in 
Islamic decoration.— Like the CFS-P 6522 structure, the 
new P4i2i2 structure can be viewed as made of inter¬ 
connected spiral chains of bonds running both clockwise 
and anti-clockwise along the z axis, with some portions of 
open space without clearly defined cavities (as, for exam¬ 
ple, in the known clathrate structures). In P4i2i2 there 
is one type of threefold spiral, three types of fourfold 
spirals and one type of fivefold spiral (whereas in CFS- 
P 6522 there is one type each of threefold, fourfold and 
sixfold spirals). The presence of fivefold spirals (Sil-Si2- 
Si2-Sil-Si3) in P4i2i2 is a rather unusual feature, which 
is unique and characteristic of this structure among oth¬ 
ers of group-14 elements previously reported or that we 
have found in our searches. There is also a large pro¬ 
portion of fivefold rings (of which there are two types) 
as well as sixfold rings. Angular distortions from the 
ideal tetrahedral angle are considerable, with the max¬ 
imum angle (126.9°) appearing along the fivefold spiral 
chains (Sil-Si3-Sil) and the minimum angle (99.7°) for 
Sil-Si 2 -Si 2 . Distances range from 2.35 to 2.39 A (cf. the 
value of 2.37 A for the diamond phase). The response of 
bonds to compression appears to be similar to diamond 
with just one bond having a somewhat larger compress¬ 
ibility. There is a large variation in the orientation of 
these bonds. 

Both CFS-P 6522 and CFS-P4i2i2 (which we could re¬ 
name CFS -6 and CFS-5, respectively, on account of the 
order of the main spirals, in each case) have rather sim¬ 
ilar energies, but the equilibrium volume of CFS-P6s22 
is larger as it is based on sixfold instead of fivefold spi¬ 
rals, which results in a more open structure, with larger 
portions of empty space. Though they have not appeared 
during our structural searches, one may legitimately won¬ 
der whether other chiral structures of the same kind but 
based upon sevenfold or even higher order spirals would 
be energetically competitive and feasible as stable phases. 


D. Band structures 

The band structures and electronic densities of states 
(EDoS) of the Pbam, P4i2i2, and P 42 /ncm structures 
of Si at zero pressure are shown in Fig. [HJ and simi¬ 
lar data for C and Ge are provided in the Supplemen¬ 
tal Material . 36 The calculated band gaps are expected to 
be underestimated, as is usual in DFT calculations using 



FIG. 7. (Color online) A slab of the P4i2i2 structure of Si 
at zero pressure, viewed along the direction of its fourfold 
rotation axis. A modicum of perspective has been introduced 
to facilitate visualization of the connectivity among sites. The 
three inequivalent sites are labelled and shown in different 
colours. 


TABLE II. The equilibrium volume per atom, Vo, bulk modu¬ 
lus, Bo, pressure derivative of the bulk modulus, B' 0 , and the 
difference in energy from the corresponding diamond phase 
at zero pressure, A E, for various structural phases of carbon. 
These calculated results were obtained with the VASP code and 
PBE density functional. (See also Supplemental Material for 
LDA and PBEsol results, as well as other phases.) 


Carbon 

Vo (A 3 ) 

Bo (GPa) 

B'o 

A E (meV) 

cd 

5.708 

435 

3.51 

0 

hd 

5.722 

433 

3.60 

25 

Pbam 

5.873 

416 

3.60 

69 

P4i2i2 

5.897 

413 

3.61 

132 

P42/ncm 

5.971 

403 

3.64 

108 

Cmca 

5.874 

413 

3.72 

128 

Cmmm 

5.870 

413 

3.67 

129 

P6 5 22 

6.232 

389 

3.51 

110 

I4i/a 

5.696 

346 

4.18 

858 

bc8 

5.616 

386 

3.84 

697 

r8 

5.658 

363 

3.98 

819 

stl2 

5.650 

395 

3.74 

886 

Ibam 

5.771 

372 

3.64 

956 

cintet 

5.498 

428 

3.72 

1115 

clat34 

6.602 

372 

3.58 

72 

clat46 

6.546 

368 

3.61 

106 


semi-local density functionals such as PBE, but the dif¬ 
ferences between the calculated band gaps are expected 
to be more accurate, by which useful estimations can be 
obtained. 

The C-Pbam, C-P4 1 2 1 2, and C-P 42 /ncm phases are 
insulators with minimum band gaps calculated to be 4.57, 
4.70, and 3.74 eV, respectively (cf. the calculated band 




















































10 



TABLE III. 

As Table [III but for silicon. 

Silicon 

Yo (A 3 ) 

Bo (GPa) 

B'o 

A E (meV) 

cd 

20.45 

88.9 

4.1 

0 

hd 

20.43 

88.6 

4.4 

11 

Pbam 

20.98 

85.0 

4.2 

29 

P4i2i2 

20.87 

85.0 

4.1 

41 

P4 2 /ncm 

21.11 

83.2 

4.4 

45 

Cmca 

20.96 

84.5 

4.2 

61 

Cmmm 

20.95 

84.8 

4.2 

61 

P6 5 22 

22.18 

80.9 

3.0 

53 

I4i/a 

18.17 

77.1 

4.4 

165 

bc8 

18.45 

83.6 

4.2 

159 

r8 

18.22 

78.8 

4.3 

160 

st 12 

18.35 

71.3 

3.1 

166 

Ibam 

18.30 

59.1 

2.8 

204 

cintet 

18.12 

96.0 

4.4 

203 

clat34 

23.56 

76.1 

4.3 

52 

clat46 

23.27 

75.8 

4.2 

67 


TABLE IV. As Table El but for 

germanium. 

Germanium 

Yo (A 3 ) 

Bo (GPa) 

B'o 

A E (meV) 

cd 

24.18 

59.2 

4.6 

0 

hd 

24.13 

59.5 

4.6 

17 

Pbam 

24.83 

56.2 

4.8 

31 

P4i2i2 

24.58 

57.5 

4.8 

39 

P4 2 /ncm 

24.82 

56.7 

4.7 

35 

Cmca 

24.81 

56.1 

4.7 

55 

Cmmm 

24.81 

56.3 

4.6 

54 

P6 5 22 

26.04 

53.2 

4.7 

33 

I4i/a 

21.74 

52.8 

5.0 

146 

bc8 

22.22 

56.5 

4.4 

137 

r8 

21.83 

53.6 

4.9 

142 

st 12 

21.73 

49.2 

4.7 

139 

Ibam 

21.92 

41.0 

5.0 

165 

cintet 

21.94 

61.7 

4.7 

173 

clat34 

27.50 

51.9 

4.8 

25 

clat46 

27.10 

51.5 

4.8 

32 


gap of C-cd of 4.13 eV). The occupied bandwidths are 
20.8 eV (C-Pbam), 19.7 eV (C-P4i2i2), and 19.9 eV (C- 
P42/ncm), which are somewhat smaller than the calcu¬ 
lated bandwidth of the diamond phase, 21.4 eV (cf. the 
experimental value of of the bandwidth of C-cd, 23.0±0.2 
eV, Ref 4^) . 

The corresponding minimum band gaps for the Si 
structures are 0.86 eV (Si-Pbam), 1.23 eV (Si-P4i2i2), 
and 0.80 eV (Si-P42/ncm), which are larger than the gap 
of 0.63 eV calculated for the diamond structure using 
PBE-DFT, which is in turn approximately half the ex¬ 
perimental value of 1.17 eV. If we add to the band gaps of 
Si-Pbam, Si-P4i2i2, and Si-P4 2 /ncm a correction equal 
to the difference between the theoretical and experimen¬ 
tal band gaps of Si-cd we obtain band gaps of about 
1.4 eV (Si-Pbam), 1.77 eV (Si-P4i2i2), and 1.34 eV (Si- 
P42/ncm). These gaps are larger than in Si-cd, and are 
considerably closer to the values around 1.5 eV that are 
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FIG. 8. (Color online) Electronic band structures and densi¬ 
ties of states of the Pbam, P4i2i2, and P42/ncm structures 
of Si at zero pressure. The valence bands are shown in red 
and the conduction bands in blue, and the band gap regions 
are indicated by parallel black lines. The zeros of energy are 
at the top of the valence band at the zone center V point. 
The panel on the right of each plot shows the total density of 
states. 


optimal for photovoltaic applications. 

In particular, the top of the valence band of Si-Pbam is 
located at the zone center T, where the (corrected) direct 
band gap is ~1.4 eV, and therefore this phase could be 
of technological interest. The bottom of the conduction 
band has a rather flat variation around T, along the T-X 
and T-Y directions (with energy differences of the order 
of 0.01 eV). Due to the small axial ratios c/a and b/a 
of the Pbam cell, both the X and Y points are close to 
the zone center in Si-Pbam, see Fig. [8] The very small 
dispersion of the conduction bands around T suggests 
that the joint density of states for electron excitation is 
large. In contrast, the band structures of both Si-P4i2i2 
and Si-P42/ncm show distinct indirect band gaps, as Si- 
cd. 

The Ge-Pbam phase has a calculated direct band gap 
at the zone center close to zero, while Ge-P4i2i2 and 
Ge-P42/ncm have calculated indirect gaps of about 1.00 
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eV and 0.30 eV, respectively. Analogous calculations for 
Ge-cd give a band gap of approximately zero (cf. the 
experimental value of 0.67 eV), thus the corrected gaps 
for the tetrahedral polymorphs are also expected to be 
larger than for the Ge-cd phase, as found in Si. 

A common feature of the band structures of the low- 
density polymorphs of group 14 phases is a weakening 
of the sp hybridization arising from the substantial de¬ 
viations of the bond angles from the perfect tetrahedral 
angle of 109.5°, and a small reduction in the occupied va¬ 
lence bandwidth. The EDoS of the Si-Pbam, Si-P4i2i2, 
and Si-P42/ncm phases are not divided into separate s 
and p parts, as found in P6s22, see Fig. 2 of Ref. |42], al¬ 
though the overlaps of the s and p parts of the EDoS in 
the Si structures are not large. The EDoS of Ge-P4i2i2 
is broken into a lower s part and an upper p part, with a 
small inner gap clearly visible at about -3.7 eV indicating 
significant weakening of the sp hybridization. The band 
structure of Ge-P42/ncm is almost broken into s and p 
parts, while for Ge-Pbam there is significant overlap in 
the energy scale between both parts, giving a gapless va¬ 
lence EDoS. 


IV. VIBRATIONAL PROPERTIES AND 
DYNAMICAL STABILITY 

The phonon dispersion relations of the Pbam, P4i2i2, 
and P42/ncm phases of Si at zero pressure, calcu¬ 
lated using a supercell method with small atomic 
displacements^ are plotted in Fig. [9j while those of C 
and Ge are given in the Supplemental Material— £ These 
plots show that the phases are all dynamically stable at 
zero pressure. The Pbam phases have distinct and rela¬ 
tively well resolved upper-frequency phonon bands which 
have rather small dispersion across the Brillouin zone, 
resulting in a characteristic high-frequency peak in their 
respective density of states (PhDoS). The phonon band 
structures and PhDoS of the three phases in Si and Ge 
phases are quite similar, differing mainly by a scaling aris¬ 
ing from the different atomic massses, which reflects the 
similarity in bonding. For the three materials, the Ph¬ 
DoS of the P42/ncm phase extends to somewhat lower 
frequencies than the Pbam phase. 

We have also simulated the Raman intensities for the 
three phases at 300 K, which are shown in Fig. [10] Si- 
Pbam has strong high-frequency Raman peaks in the 
range 450-510 cm -1 , and a number of rather weak Raman 
modes below that frequency range. Its strongest mode 
lies almost directly on top of the calculated single Raman 
mode of Si-cd (505 cm -1 , cf. the experimental value of 
520 cm" 1 , Ref.0). This fact is related to the structural 
similarity of Pbam and diamond discussed previously. In 
contrast, the strongest Raman peak of Si-P4i2i2 lies at 
a frequency of 480 cm -1 , somewhat below the diamond 
mode. Si-P4i2i2 shows also significant Raman activity 
just below 200 cm -1 and weak Raman activity around 
330 cm -1 , as can be seen in Fig. [T0l 



FIG. 9. (Color online) Calculated phonon dispersion relations 
of the Pbam, P4i2i2, and P42/ncm structures of Si at zero 
pressure. The panel on the right of each plot shows the total 
phonon density of states. 

These simulated spectra can be used for the eventual 
identification of the phases in Raman experiments, al¬ 
though insufficient experimental resolution, broadened 
peaks, and mixture of phases with overlapping peaks (all 
of them rather normal situations in nanoindentation or 
conventional compression experiments) difficult the task 
of comparison with experimental Raman results. Also, 
some caution has to be exerted as changing temperature 
can result in a significant change in the relative inten¬ 
sity of the calculated diffraction peaks. For example, 
we find significant differences in the relative intensities 
of the modes reported for Si-P42/ncm by Zhao et al.j& 
probably due to the temperature used in simulating the 
Raman spectra. (Further to this, Zhao et al. used the 
LDA while the results in Fig. [TO] were obtained using the 
PBE functional.) 

Zhao et also made the interesting proposal that 
P42/ncm could explain the experimentally unknown 
structure of a Si-XIII phase, supported by a compari¬ 
son of calculated and experimental Raman vibrational 
data.—! i L 21 i 22 Metastable Si-XIII is relevant to thermal 
processing of silicon wafers, which induces large and po- 
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Frequency (cm *) 

FIG. 10. (Color online) Calculated Raman spectra of the 
P42/ncm, P4i2i2, and Pbam structures of Si at zero pressure 
and 300 K, calculated within the PBE. The vertical dotted 
line indicates the calculated frequency of the single intense 
Raman mode of the Si-cd phase. 


tentially damaging stresses. It has been reported in sev¬ 
eral source a 6 i 7 i 20 ~ — as being observed only in coexistence 
with other Si phases (bc8, r8, cd, and hd), which makes 
it very difficult to determine its structure. The experi¬ 
mental situation concerning Si-XIII has been examined 
in detail in a recent publication by Ruffell et al£ who 
have carefully reassessed previous experimental data on 
this phase^flr— These authors conclude the existence 
of three signature peaks which are characteristics of the 
Raman spectrum of Si-XIII, at frequencies of 202, 333, 
and 478 cm” 1 , which agree quite well with those calcu¬ 
lated by us for the P4i2i2 phase (cf. the data in Fig. 
[TO] for Si-P4i2i2 with those in Fig. 6 of Ref. 0 for the 
Si-XIII phase). In contrast, the Raman spectrum of the 
P4 2 /ncm structure, previously suggested for Si-XIII by 
Zhao et al.j& lacks the characteristic peak at 333 cm” 1 , 
whereas their calculated peaks at ~110 and ~380 cm” 1 
are absent in the experimental data for this phase. We 
suggest thus that P4i2i2 is a more likely candidate for 
the structure of the Si-XIII phase. 


V. SUMMARY AND CONCLUSIONS 

Structure searches for group 14 elements using AIRSS 
have led to the discovery of a number of low-energy 
structures, including low-density sp 3 bonded structures 
of Pbam and P4i2i2 symmetries. From the structures 
that we have studied in C, Si, and Ge, only the cd and hd 
structures have significantly lower energies than Pbam. 
Another structure that readily appeared in our searches 
is the P4 2 /ncm (or T12) structure previously reported 
by Zhao et alM* for which we provide a detailed interpre¬ 
tation in terms of tilting and rebonding of diamond-like 
tetrahedra, that complements and extends those authors’ 
description. The P4 2 /ncm structure is the simplest pos¬ 
sible that can be built by stacking slabs of tilted tetra¬ 
hedra while preserving to a large degree a highly regular 
fourfold coordination for the sites. The Pbam, P4i2i2, 
and P4 2 /ncm structures are dynamically stable and have 
lower energies than a number of metastable polymorphs 
of group 14 elements that have been synthesized. It is 
conceivable that one or more of these phases could be 
formed as metastable phases in indentation or diamond 
anvil cell experiments, perhaps in combination with heat 
treatment, and in fact P4 2 /ncm has been proposed previ¬ 
ously as the structure of known yet experimentally unre¬ 
solved phases of Si and Ge^ We note that an equivalent 
P4 2 /ncm structure has recently been suggested to play 
an important role in the homogeneous crystallization of 
water 

Our PBE-DFT calculations show semiconducting be¬ 
havior for the Pbam and P4i2i2 structures at zero pres¬ 
sure, except for Ge-Pbam that is calculated to have a 
band gap close to zero at low pressures. The band gaps 
for a particular structure decrease from C to Ge, as ex¬ 
pected. The direct band gap of Si-Pbam at the zone 
center is estimated to have a value of 1.4 eV (after cor¬ 
recting for the well-known underestimation of the PBE- 
DFT gap) so that it might be suitable for applications in 
photovoltaics. 

We have found several other structures in our searches. 
These include a low-energy structure of J4i/a symmetry 
that is denser than the cd structure in Si and has a very 
similar enthalpy to that of the well-known metastable 
polymorphs r8/bc8, and which has been suggested to 
play a role in the decompression kinetics from the high 
pressure /3-Sn phased We also found a low energy and 
low density structure of Cmca symmetry which, in car¬ 
bon, is very close in energy to both the Pmmn struc¬ 
ture of the so-called P phased and the Cmmm structure 
or Cco-C8 phased with which it shares similar struc¬ 
tural features. These phases have been proposed for su- 
perhard carbon allotropes experimentally obtained after 
cold compression of carbon nanotubes^i and their close¬ 
ness suggests that such experimental allotropes could 
consist of a mixture of energetically and structurally re¬ 
lated forms. 

Among the proposed structures, the new Pbam poly¬ 
morph looks particularly promising (low energy, dynam- 
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ical stability and direct band gap in Si). Although 
the study of transformation mechanisms between phases 
and possible synthesis routes is beyond the scope of the 
present work, its structural features, based on sheared 
diamond-like slabs, suggests that it might be synthesized 
by controlled application of uniaxial stress to the dia¬ 
mond phases. (Note added: When this paper was about 
to be submitted we learnt of very recent findings for 
Pbam in C by Baburin et air 7 Our results for C-Pbam 
agree with those in Ref. while here we also provide 
results for Si and Ge4£) 

P4i2i2, a spiral structure with tetragonal symme¬ 
try consisting of a packing of pentagonal helices with 
a large proportion of fivefold rings, is also very inter¬ 
esting and shows several unusual features. It belongs 
to the same class of chiral frameworks as a previously 
reported P 6522 structure^ and its projection along the 
four-fold axis yields the so-called Cairo two-dimensional 
pentagonal tiling. In a very recent publication^ that ap¬ 
peared when this paper was about to be submitted, a 
new two-dimensional carbon allotrope with 24 atoms per 
cell and P42im symmetry was proposed, with exactly 
the same description in terms of the Cairo tiling as the 
three-dimensional P4i2i2 structure proposed here. The 
appearance of the same unusual (yet beautiful) topologi¬ 
cal motif in both two- and three-dimensional metastable 
forms of carbon, obtained independently, is quite inter¬ 
esting and begs for a closer examination of the role of pen¬ 


tagonal rings on the stability of novel and exotic carbon 
allotropes. The P4i2i2 structure is energetically quite 
competitive in Si and Ge, and in C it has a similar en¬ 
thalpy to the structures currently proposed for allotropes 
obtained from compressed nanotubes j 27 ' 30 i 31 yet with a 
very different topology. We find that the P4i2i2 struc¬ 
ture shows the best compatibility with the available ex¬ 
perimental Raman data for the unknown phase Si-XIII. 

We hope that these results will entice experimental 
searches for new phases of these materials. 
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